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ABSTRACT 

Molecular clouds are expected to emit non-thermal radiation due to cosmic ray inter- 
actions in the dense magnetized gas. Such emission is amplified if a cloud is located 
close to an accelerator of cosmic rays and if energetic particles can leave the acceler- 
ator site and diffusively reach the cloud. We consider here the situation in which a 
molecular cloud is located in the proximity of a supernova remnant which is efficiently 
accelerating cosmic rays and gradually releasing them in the interstellar medium. We 
calculate the multiwavelength spectrum from radio to gamma rays which is emerging 
from the cloud as the result of cosmic ray interactions. The total energy output is 
dominated by the gamma ray emission, which can exceed the emission in other bands 
by an order of magnitude or more. This suggests that some of the unidentified TeV 
sources detected so far, with no obvious or very weak counterparts in other wave- 
lengths, might be in fact associated with clouds illuminated by cosmic rays coming 
from a nearby source. Moreover, under certain conditions, the gamma ray spectrum 
exhibit a concave shape, being steep at low energies and hard at high energies. This 
fact might have important implications for the studies of the spectral compatibility of 
GeV and TeV gamma ray sources. 
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1 INTRODUCTION 

Galactic Cosmic Rays (CRs) are believed to be ac- 
celerated via first order Fermi mechanism at the ex- 
panding shock waves of supernova remnants (SNRs) 
(for a review see iDrurv 19831 ; iBlandford and Eichler 19871 ; 
I Jones and Ellison 199ll ). Two of the strongest points sup- 
porting this scenario are: i) the fact that SNRs alone 
are capable to maintain the galactic CR flux at the ob- 
served level, provided that a fraction of about 10% of 
their explosion energy i s somehow converted into CRs 
l|Baade and Zwickv 19341 ; iHillas 20051 ). and ii) the fact that 
som e shell-type SNRs hav e been detected at TeV energies 
(see I Aharonian et al. 20081 . for a review), as expect ed if they 
indeed are the sources of CRs l|Drurv et al. 19941 ). In this 
context, TeV gamma rays originate from the decay of neu- 
tral pions produced in interactions between the accelerated 
CRs and the interstellar gas swept up by the SNR shock. 

However, the detection of SNRs in gamma rays, though 
encouraging, cannot provide by itself the final proof that 
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CRs are indeed accelerated at SNR shocks. This is be- 
cause competing leptonic processes, namely inverse Comp- 
ton scattering from accelerated electrons, can also explain 
the observed TeV emission, provided that the magnetic field 
does not significantly exceed ~ 10 [iG (e.g. (author?) 
I Aharonian el al. 2006a! ) . Evidence for strong ~ 100 /iG mag- 
netic field, and thus indirect support to the hadronic sce- 
nario for the gamma ray emission, comes from the ob- 
servation of th in X-ray synchrot r on filaments surround- 
ing some SNRs feamba et al. 20031 ; I Vink and Laming 20031 ; 
IVolk et al. 20051 ) and of the ra pid variability time 
scale of the synchrotron X- rays (|Uchivama et al. 20071 ; 
lUchivama and Aharonian 20081 ). On the other hand, the 
absence of thermal X-ray emission from the TeV SNRs 
RXJ1713. 7-3946 and RXJ0852.0-4622 has been interpreted 
by Katz and Waxman (2008) as an argument against 
the hadronic nature of the gamma-ray emission (and 
thus supporting the leptonic scenario). Their argument 
can be summarized as follows: in order to fit the TeV 
data with an hadronic model, the gas density can- 
not be too low. The high density would then imply 
a strong X-ray thermal Bremsstrahlung emission which 
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is not observed. However, this argument is not robust 
given the uncertainties related to the heating of elec- 
trons at shocks, and it has been st r ongly questioned 
dTanaka et al 20081; iMorlino et al. 20081 ; IVolk et al. 20081 ; 
iGabici 20081 ; iDrurv et al. 20081 ). Moreover, it has been 
shown that the downstream proton temperature can be ef- 
fectively reduced if the shock is effectively accelerating cos- 
mic rays. This would also suppress th e electron temper- 
ature and the the rmal X-ray emission (|Drurv et al. 20081 ; 
iHelder et al. 20081) . A conclusive proof of the hadronic na- 
ture of the gamma-ray emission will possibly come from the 
detection of neutrinos from SNRs, which are expected to be 
produced during the same hadronic in teractions responsi- 
ble for the production of gamma rays l)Kappes et al. 20071 ; 
ICostantini and Vissani 20051 ; IGabici and Aharonian 2007h . 
However, since the detection of neutrinos appears challeng- 
ing even for km 3 -scale neutrino telescopes, the search for 
evidence of CR proton acceleration coming from gamma ray 
or multi wavelength observations is mandatory. 

The presence of a massive molecular cloud close to a 
SNR can provide a dense target for CR hadronic interac- 
tions and thus enhance the expected gamma-ray emission. A 
correlation between SNRs and mo lecular clouds is e xpected, 
especially in star forming regions ()Montmerle 19791 ). In gen- 
eral, a spatial correlation between TeV gamma rays and 
dense gas, though not conclusive, would favor an hadronic 
scenario, where gamma rays are the result of the inter- 
actions of CR protons in the dense gas. Such correlation 
has been obser ved from a number of S N R/molecular clouds 
associations ( Aharonia n et al. 2008b]; lAlbert et al. 20071 : 
Ah aroman et al. 2008al ). If the molecular cloud is overtaken 
by the SNR shock, the enhanced gamma-ray emission is ex- 
pecte d to be cospatial with the SNR shell, or with a portion 
of it (I Aharonian et al. 19941 ). If the cloud is located at some 
distance from the SNR, it can still be illuminated by CRs 
that escape from the SNR and produce gamma rays ther e 
(I Aharonian and Atovan 19961 ; IGabici and Aharonian 20071 ). 
For this scenario, it has been shown that, for typical SNR 
parameters and for a distance D = 1 kpc, a molecular cloud 
of mass 10 4 M Q can emit TeV gamma rays at a detectable 
level if it is located within few hundred parsecs from the SNR 
ijGabici and Aharonian 20071 ). This implies that the angu- 
lar displacement between the SNR shell and the gamma ray 
emission is of the order of d « 6°(D/lKpc)- 1 (d c; /100pc), 
where d c i is the distance between the SNR and the cloud. 
This translates in the fact that sometimes the association 
between SNRs and molecular clouds can be not so obvi- 
ous, given that the separation between the two objects can 
be even larger than the detector's field of view. Following 
this rationale, Gabici and Aharonian (2007) proposed that 
some of the unidentifie d TeV sour ces detected by HEGRA 
|Aharonian et al. 20021). H.E.S.S. l|Aharonian et al. 2006b! ) 



and MILAGRO (|Abdo et al. 20071 ) might in fact be such 
clouds illuminated by a nearby SNR. 

In this paper, we calculate the expected non-thermal 
emission, from radio to multi- TeV photons, from a molecu- 
lar cloud illuminated by CRs coming from a nearby SNR. 
We generalize the model presented by Gabici and Aharo- 
nian (2007), which was limited to the hadronic TeV photons 
only, to include the generation of secondary electrons in the 
cloud and the related synchrotron and Bremsstrahlung emis- 
sion. We found that the total radiation energy output from 



a cloud is dominated by the gamma ray emission, which 
can exceed the emission in other bands by an order of mag- 
nitude or more. This suggests that some of the unidenti- 
fied TeV sources detected so far, with no obvious or very 
weak counterparts in other wavelengths (the so called "dark 
sources"), might be in fact associated with massive molecu- 
lar clouds illuminated by CRs. Moreover, under certain con- 
ditions, the gamma-ray spectrum from the cloud exhibit a 
concave shape, being steep at low (~ GeV) energies and hard 
at high (~ TeV) energies. This fact might have important 
implications for the studies of the spectral compatibility of 
GeV and TeV gamma ray sources. 



2 THE MODEL 

Consider a supernova of total energy 10 51 i?5ierg explod- 
ing in a medium of density n. The initial shock velocity is 
10 9 ug cm/s and remains roughly constant until the mass of 
the swept up material equals the mass of the ejecta. This 
happens at a time tsedov ~ 200[£ , 5i/(n ug)] 1//3 yr, when the 
shock radius is tti 2.1[E 51 /(n w|)] 1/3 pc. Then the SNR en- 
ters the Sedov phase and the shock radius and velocity scales 
with time as R s h oc t 2 ^ 5 and u s h cx i~ 3//5 . 

The spectrum of particles accelerated at t he SNR shoc k 
is determined by the transport equation fe.g.. IDrurv 19831 ): 



f-V-/ + uV/-^g=0. 



(1) 



where D — D(p) is the momentum dependent diffusion co- 
efficient and u the flow velocity. For a strong shock with 
compression factor r s = 4, the test particle theory predicts 
an u niversal shap e for the CR spectrum at the shock fo (p) oc 
p~ A (|Drurv 19831 ). If the shock is an efficient accelerator 
(as SNR shocks are believed to be), the CR pressure mod- 
ifies the flow structure, making the shock more compress- 
ible and the spectrum of the ac celerated particles harde r, 
fo(p) oc p~ a with 3.5 < a 4 l|Malkov and Drurv 20o"ll ). 
Detailed calculations compared with multiwavelength ob- 
serv ations of SNRs suggest the values r s ~ 7 and a w 
3.7 |Ellison et al. 20071 ; iBerezhko and Volk 20061 ). which we 
adopt in the following. 

The maximum momentum of the accelerated particles 
is determined by a confinement condition, namely that the 
diffusion length Id of the particles cannot exceed the char- 
acteristic size of the system R sn '- 



, D(p max ) ^ n 

'd — ^ fish- 

Ush 



(2) 



The maximum possible energies are achieved when the ac- 
celeration proceeds in the Bohm diffusion limit, D oc p/B s h, 
with B s h the magnetic field strength at the shock. In 
this case the maximum momentum decreases with time as 
Pmaxit) oc Bght' 1 ' 5 . In fact, the drop of Pmax is even faster, 
given that the magnetic field is also expected to decrease 
with time. This implies that at any time, particles with mo- 
mentum above p m ax(t) quickly escape the remnant, gener- 
ating a cutoff in the spectrum. The spectrum of the run- 
away particles can be approx imated as a 5-function (see 
IPtuskin and Zirakashvili 20051 ): 

q ESC {p,t) = -8(p - Pmax) 
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where the integration has to be performed where the inte- 
grand is negative. Thus, to calculate the flux of the run- 
away particles one has to know: (i) the CR particle distri- 
bution function at p max at any location in the SNR, (ii) the 
flow velocity both inside the shock and outside it, where the 
CR precursor forms and (Hi) how the maximum momentum 
varies during the SNR evolution. Ptuskin and Zirakashvili 
(2005) showed that it is straigthforward to derive (i) and (ii) 
using an approximate (but still reasonably accurate) linear 
velocity profile inside the SNR : 

u s h{t) 



U= I — 



Rsh{t) 



R 



(4) 



and assuming that the CR pressure at the shock P^" is a 
fraction £cr of the incoming ram pressure gu 2 3h and that 

foiPrnax) OC Pg* . 

The determination of p ma x and its evolution with time 
requires the knowledge of the diffusion coefficient (see Eq. 
[2| , which is in turn determined by the level of magnetic tur- 
bulence generated by the accelerated particles themselves. 
This makes the problem nonlinear and very difficult to be 
solved. The value of p ma x depends on a few crucial but 
poorly known aspects of the problem, including the na- 
ture of CR-driven instability operati ng in the shock pre- 
cursor and the level of wave damping fBell and Lucek 2001: 
Bell 2004 IPtuskin and Zirakashvili 20051 ; iBlasi et al. 20071 ; 
Vladimirov et al. 20061 ). Because of these uncertainties, 



we adopt here a phenomenological approach, namely we 
parametrize the maximum momentum as p max (t) oc t~ s . 
We further assume p m ax ~ 5 PeV and ~ f GeV at the 
early (t = 200 yr) and late (t — 5 x f0 4 yr) epochs of 
the Sedov phase respectively. This requires S ~ 2.48. Re- 
markably, if the maximum momentum is a power law func- 
tion of time, the spectrum of the escaping particles inte- 
grated over t he whole Sedov phase is also a p ower law of the 
form oc p~ 4 (|Ptuskin and Zirakashvili 20 05). which is close 
(sligthly harder) to what needed to fit the CR data below 
the knee (jBerezinskii et al. 1990l ). 

Following Ptuskin and Zirakashvili (2003), the approxi- 
mate spectrum of the CRs inside the SNR /j n (R, p, t) can be 
obtained from Eq. [T] by dropping the diffusion term, while 
the spectrum of the runaway CRs at a given distance R from 
the SNR and at a given time t can be obtained by solving 
the diffusion equation: 



(R,p,t) = D IS m(p)V fout(R,P,t) + q eS c(p,t)5(R).(5) 



dfout 

~dT 

The diffusion coefficient Dism(p) describes the propaga- 
tion of CRs in the galactic disk. The available CR data re- 
quire a power-law energy dependence, Dism(E) cx E~ s , 
with D IS m ~ 10 28 cm 2 /s at E « fO GeV and s w 
0.3 0.7 (|Berezinskii et al. 1990h . The constraints on the 
diffusion coefficient are obtained from the comparison be- 
tween diffusion models and CR data and have to be con- 
sidered as average galactic values. However, the conditions 
might be rather different in regions close to CR sources, 
in particular due to the presence of strong gradients in 
the CR distribution, which may enhance the generation 
of plasma waves and thus suppress the diffusion coefficient 
|Wentzel 19741 ; IPtuskin et al. 2 008). The change in s within 



the allowed range or the choice of a different normalization 
for Dism does not alter qualitatively the results, the main 
effect being that the characteristic time scales of the problem 
change proportional to 1/Dism- 

Remarkably, if p m ax (t) scales as a power law of time, Eq. 
[5] can be solved analytically and the distribution function of 
escaping cosmic rays at any given distance R from the SNR 
and at any given time t reads, for energies E ^ cx p max (t): 



f out {t,R,E) = 
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where Esn is the total supernova explosion energy, rj is the 
fraction of such energy converted into CRs and Em ax and 
Emin are the maximum and minimum energies of CRs ac- 
celerated during the Sedov phase. The diffusion distance for 
a CR of energy E is: 



R d (E) = y/AD(E)(t- x (E)) 



where 



X(E) = ts e do 



•(ft-) 

V &MAX I 



-1/6 



(7) 



(8) 



> Em ax , 

represents the time after the supernova explosion at which 
CRs with energy E are released in the interstellar medium. 
The solution derived by Ptuskin and Zirakashvili (2005) for 
the total CR spectrum injected by a SNR in the interstellar 
medium during the whole Sedov phase can be easily derived 
by integrating Eq.[6]over space. Finally, it has to be noticed 
that the total CR spectrum at a given time and at a given 
distance from the SNR is the sum of two contributions: i) 
a time dependent contribution from CRs coming from the 
SNR, whose spectrum is described by Eq.[6]and ii) a steady 
contribution from the galactic CR background. 

Following the procedure described above, it is possible 
to evaluate, for any given time, the CR spectrum in prox- 
imity of a molecular cloud located at a given distance from 
the SNR. If the diffusion coefficient inside the cloud is not 
significantly smaller than the Galactic one, CR freely pene- 
trate the cloud and the CR spectrum inside the cloud is not 
affected by propagation effects. Conversely, if the diffusion 
coefficient is significantly (more than one order of magni- 
tude) reduced, low energy CRs are excluded from the cloud 
and a low energy cutoff appears in the CR spectrum, at 
an energy that depends on t he value of the diffusion coef- 
ficient (see lGabici et al. 2 007. for details). Here, we assume 
free penetration of CRs and we refer the reader to Gabici 
et al. (2007) (and references therein) for a detailed discus- 
sion on CR exclusions from molecular clouds. We do not 
consider here any contribution from CR electrons coming 
from the SNR, since they do not escape the remnant due 
to diffusive confinement (for low energy electrons) and se- 
vere synchrotron losses in the strong magnetic field (for high 
energy electrons). 

CR protons propagating inside a molecular cloud pro- 
duce secondary electron-positron pairs during inelastic in- 
teractions in the dense intercloud medium. We calculate the 
spectrum of the injected secondary electrons Q e (t,E) by 
using the parameterization given by Kelner et al. (2006) 
and we follow the time evolution of the electron distribution 
function f E (t,E) by solving the equation: 



dfe(t,E) 

dt 



d_ 

dE 
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+ Qe(t,E) 



Mt,E) 



(9) 
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where (dE/dt) B = E/rf osa is the energy loss rate for elec- 
trons, rf oss the energy loss time, and tI sc the diffusive escape 
time from the cloud. All these time scales will be discussed 
and estimated in section 4. Once the proton and electron CR 
spectra have been derived, the non-thermal radiation from 
the molecular cloud can be calculated. 



3 COSMIC RAY SPECTRUM AT THE CLOUD 
LOCATION 

We consider here a molecular cloud located at a given dis- 
tance d c i from a SNR and we calculate the CR spectrum at 
the cloud location. As mentioned above, the spectrum is the 
sum of two disctinct components: i) the CRs coming from 
the nearby SNR, described by Eq. |SJ and ii) the galactic 
CR background, which results from the superposition of all 
the CR sources in the Galaxy. While the latter contribution 
is constant in time, the first one changes, since the flux of 
CRs escaping from the SNR evolves in time as described in 
SecH 

Fig- Q] shows the spectrum of CRs at the location of the 
molecular cloud. The Galactic CR background is plotted as 
a thin dot-dashed line labeled as CR sea, while the spec- 
trum of the CRs coming from the SNR is plotted as a thin 
line for different times after the supernova explosion: 500 yr 
(solid), 2000 yr (dotted), 8000 yr (short - dashed) and 32000 
yr (long-dashed). Thick lines represent the sum of the two 
contributions. The distance between the SNR and the molec- 
ular cloud is 50 pc (left panel) and 100 pc (right panel). We 
assume a total supernova explosion energy of 10 51 ergs and 
a high CR acceleration efficiency at the SNR shock equal 
to 77 = 30%. The normalization of the CR spectrum at the 
cloud scales linearly with these two quantities Q. The diffuse 
Galactic cosmic ray spectrum is assum ed to be ident ical to 
the one observed near the Sun (see e.g. iDermer 1 986): 

Jcr(E) = 2.2 (j^y) cm-V^eV^sr- 1 (10) 

while the diffusion coefficient, needed to evaluate Eq. [6] is 
taken equal to: 

D w = 1028 ( idb) °' 5 cm2/a ' 

compatible with CR propagation models (e.g. 
iBerezinskii et al. 1990h . 

The evolution with time of the CR spectrum at the po- 
sition of the molecular cloud can be understood as follows. 
According to the model described in Sec. [2] CRs with differ- 
ent energies leave the SNR at different times. The highest 
energy (~ PeV) CRs leave the SNR first, while CRs with 
lower and lower energy are released at later times. More- 
over, higher energy CRs diffuse faster, thus the spectrum of 
CRs at the cloud exhibit a sharp low energy cutoff at an en- 
ergy Eiow> which moves to lower and lower energies as time 
passes. The position of the cutoff represents the energy of 
the least energetic particles that had enough time to reach 
the cloud. 



From Fig. [T] it is clear that the influence of the pres- 
ence of a nearby SNR close to the cloud is reflected in the 
CR spectrum at the cloud position, but this influence de- 
pends on many parameters, such as the distance between 
the SNR and the cloud, the time since the supernova explo- 
sion, and the CR energy. The CRs coming from the SNR 
dominate the total CR spectrum at high energy, while at 
lower (~ GeV) energies the galactic CR background is al- 
ways the dominant component, unless the molecular cloud is 
located at distances significantly smaller than ~ 50 pc from 
the SNR. However, such small distances are comparable to 
the size of the SNR itself and thus in this case an inter- 
action between t he SNR shock and the molecular cloud is 
expected (see e.g. lAharonian et al. 19941 ). The investigation 
of this scenario goes beyond the scope of this paper and thus 
we limit ourselves to the case in which the distance between 
the SNR and the cloud is significantly bigger than the size 
of both objects. 

Aharonian and Atoyan (1996) also evaluated the CR 
spectrum in the vicinity of a CR accelerator by using an 
approach similar to the one developed here. They made no 
specific assumption about the nature of the accelerators and 
solved the CR transport equation by assuming that a power 
law spectrum of CRs is injected in the interstellar medium. 
They considered both the case of a continuous injection of 
particles during the whole lifetime of the accelerator and the 
case of an impulsive source that releases all the CRs at the 
same time. The approach we adopt here is different, because 
it is specific for a given class of sources, namely SNRs. In this 
specific case, particles having different energies are released 
at different times in the interstellar medium. 



4 RELEVANT TIME SCALES FOR COSMIC 
RAY PROPAGATION INSIDE A 
MOLECULAR CLOUD 

Molecular clouds are characterized by a wide range of 
masses, going from « 10 M© to 10 s M® or even more and 
have typical sizes ranging from few to few tens of par- 
sees. The typical density of a cloud is of about few hun- 
dred atoms per cubic centimeter, but much higher densi- 
ties can be found in less massive and smaller (parsec or 
sub-parsec s cale) molecular cloud co res, dark clouds or Bok 
globules fsee lStahler and Palla 20051 . for a review). The typ- 
ical magnetic fi eld of the intercloud medium is m 10 fiG 
l|Shu et al. 19871 ). and it scales roughly as the square root of 
the gas density, thus reaching the m G level in the d ensest 
regions with density 10 5 -i- 10 6 cm -3 (|Crutcher 19991 ). 

The propagation of high energy CRs inside molecular 
clouds has been studied in Gabici et al (2007) (see also Do- 
gel' and Sharov, 1990), where an extensive discussion can be 
found. Here we summarize the most relevant aspects. Once 
the CRs from the SNR reach the molecular cloud, they dif- 
fusively penetrate with typical time scale: 

R 2 

Tdiff ~ 6D(E,B) (12) 



1 This is true only as a first approxiamation for the parameter 
77, since the estimated flux of particles escaping the SNR comes 
from a nonlinear theory of shock acceleration 



where R c i is the cloud radius and D is the diffusion coef- 
ficient which we assume here to depend on energy and on 
magnetic field as: 
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Figure 1. Spectrum of CRs at the location of the molecular cloud. The cloud is located at 50 pc (left) and 100 pc (right panel) from 
a SNR. The thin dot-dashed line shows the Galactic CR spectrum, while the thin solid (curve 1), dotted (2), short-dashed (3) and 
long-dashed (4) lines represent the spectrum of CRs coming from the SNR for 500, 2000, 8000 and 32000 years after the supernova 
explosion, respectively. The thick lines show the total CR spectrum at the cloud location. 
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Here \ is a factor that takes into account deviations from 
the average Galactic di ffusion coefficient described by Eq. llll 
(see iGabici et al. 20071 ) and 3 fiG is the average magnetic 
field in the Galactic disk. At very high energies, when the 
Larmor radius of the particle becomes comparable with or 
even larger than the size of the cloud, particles propagate 
almost rectilinearly, and the propagation time reduces to 
the crossing time r cro3S = R c i /c. In the following we assume 
the characteristic propagation time for a CR in a molecular 
cloud to be: 



Tdiff +t c 



(14) 



which is a rough approximation which still describes with 
sufficient accuracy the two different regimes of propagation. 

CRs can freely penetrate the molecular cloud if the dif- 
fusion time is shorter than the energy loss time which, for 
CR protons with energy above ~ lGeV, is dominated by in- 
elastic p roton-proton interactio ns in the dense gas and reads 
(see e.g. iBerezinskii et al. 1990l ): 



6 x 10' 



\ 100 cm~ 3 ) 



yr, 



(15) 



where n gas is the gas density, c is the speed of light and 
a pp and k are the cross section and inelasticity of the pro- 
cess, respectively. For nonrelativistic protons (energies below 
1 G eV) ionization losses b ecome relevant, with time scale 
(e.g. IBerezinskii et al. 1990l ): 



2.8 x 10 



7 / "-gas \ 
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E k 



yr. 



(16) 



Here E k is the proton kinetic energy and m v is the proton 
mass. The total energy loss time for CR protons due to both 
processes is given by: 



+ 



(17) 



Secondary electron-positron pairs are produced inside 
the cloud during inelastic collisions between CR protons 
in the dense gas. The typical diffusion time for such elec- 
trons is given by Eq. 1121 while the relevant channels for 
electron energy losses are ionization losses, Bremsstrahlung 
emission a nd synchrotron emission with ch aracteristic time 
scales (e.g. iGinzburg and Svrovatskii 19641 ): 
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(18) 
(19) 
(20) 



respectively. Here, 7 is the electron Lorenz factor and E the 
total electron energy. The ionization losses are computed for 
ultra-relativistic electrons. The total energy loss time for CR 
electrons due to the three processes listed above is given by: 



7~loss 



1 



+ 



1 



1 



TBre 



(21) 



In writing the equation above, we neglected the possible role 
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Figure 2. LEFT PANEL: relevant time scales for CR propagation inside a molecular cloud with mass M c i = 10 5 Mq, radius i? c ; = 
20 pc, and magnetic field 20 fiG. Assuming a flat density profile gives a density of ~ 120 cm -3 . The dashed line represents the 
CR propagation time scale over a distance R c i. The dotted line represents the energy loss time for CR protons (ionization losses are 
relevant below 1 GeV, inelastic proton-proton interaction at higher energies) while the solid line refers to the energy loss time for CR 
electrons, including ionization losses, Bremsstrahlung losses and synchrotron losses, which dominates at low, intermediate and high 
energies respectively. MIDDLE PANEL: same as left panel, but for a cloud with mass M c i = 1O 3 M0, radius R c i = 0.5 pc, and 
magnetic field 300 /xG. The corresponding density is ~ 1.6 X 10 5 cm — 3 . RIGHT PANEL: same as the left panel, but for the specific 
case of the SgrB2 cloud with mass 1.9 X 10 6 Mq, radius 5.5 pc and magnetic field 1 mG. This implies an average density of 1.1 X 10 5 
cm -3 (see text for an explanation of the choice of parameters adopted). 



of inverse Compton losses off soft background photons. Such 
losses are expected to be relevant when the energy density 
in the radiation field u) ra d is greater than the magnetic field 
energy density, fn the following we will consider a cloud with 
radius R c i = 20 pc and magnetic field B — 20/iG, and in this 
case the condition reads: u) ra d > 10 (5/20 fiG) 2 eV/cm 3 . 
Such condition can be realized inside molecular clouds only 
close to star forming regions. In particular, if a star forming 
region or a stellar OB association with total photon output 
equal to 4 x IQ^L/Lq erg/s is located within a molecu- 
lar cloud, inverse Compton losses will dominate over syn- 
chrotron losses around the star forming region up to a dis- 
tance of R « 8 X 10~ 3 (L/L Q ) 1/2 (B/20 ^G)" 1 pc, which be- 
comes comparable with or larger than the molecular cloud 
radius R c[ when i ) 6 x 10 6 (B/20 fiG) 2 (R ci /20 pc) 2 L©. 
Since this is a quite high luminosity, only a small fraction of 
molecular clouds are expected to host such luminous OB as- 
sociations (Williams and McKee 1997). We can thus safely 
neglect the role of inverse Compton losses in most of the 
situations. Inverse Compton losses can certainly play a role 
if one considers small regions within the cloud which sur- 
rounds star forming regions. Such regions constitute a small 
fraction of the cloud volume and thus are not expected to 
affect significantly the non thermal emission from the whole 
cloud. 

In Fig. [2] (left panel) the characteristic time scales listed 
above have been plotted as a function of particle energy 
for a giant molecular cloud with total mass M cl = 1O 5 M 
and radius R c i = 20 pc. Assuming a flat density profile 
the density is n gas ~ 120 cm -3 . The magnetic field is as- 
sumed to be B c i = 20 /iG. The dotted line refers to proton 
energy losses, which are dominated by ionization losses at 
energies below ~ 1 GeV and by inelastic proton-proton in- 
teractions at higher energies. The solid line represents the 
electron energy loss time. The three different power law be- 



haviors reflect the dominance of ionization, Bremsstrahlung 
and synchrotron losses at low, intermediate and high en- 
ergies, respectively. Finally, the dashed line represents the 
propagation time over a distance equal to the cloud radius. 
The propagation time has been evaluated by using the dif- 
fusion coefficient from Eq. [T3] with \ = 1 ( no suppression 
with respect to the average Galactic value). 

For proton energies above the threshold for pion pro- 
duction (E t h ~ 280 MeV), the propagation time is always 
shorter than the energy loss time. This means that CR pro- 
tons which produce gamma rays and secondary electrons can 
freely penetrate the cloud and their flux is not attenuated 
due to energy losses. The propagation time for CR electrons 
is also shorter than the energy loss time for particle energies 
between E ~ 100 MeV and a few hundreds of TeV. This 
implies that, within this energy range, the secondary elec- 
trons produced inside the cloud quickly escape, and have 
little effect on the non-thermal emission from the cloud. On 
the other hand, extremely energetic electrons with energies 
above a few hundreds TeVs radiate all their energy in form 
of synchrotron photons before leaving the cloud. In a typical 
magnetic field of a few tens of microGauss, these electrons 
emit synchrotron photons with energy: 

E s1J n *** 1 I — — ^— I ( VkeV. (22) 

y \ 10 /j,G J V 100 TeV / V ; 

Thus, the most relevant contribution from secondary elec- 
trons to the cloud non thermal emission falls in the hard 
X-ray band. 

The middle panel of Fig. [5] shows the time scales for 
a compact cloud with radius R c i = 0.5 pc and mass M c i = 
1O 3 M0, which implies, in case of a flat density profile, a den- 
sity of ~ 1.6x 10 5 cm -3 . A strong magnetic field of 300 fj,G is 
assumed, as sug gested by observations of very dense clouds 
|Crutcher 19991 ) . These parameters are typical for molecular 
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cloud cores or for dark clouds (|Stahler and Palla 20051 ). The 
left and middle panels in Fig.[2]looks very similar, except for 
a scaling factor. This implies that the same conclusions can 
be drawn as for the case of a giant molecular cloud, and thus 
a similar behavior is expected from giant molecular clouds 
and compact clouds, the only difference being that all the 
time scales are much shorter in the latter case, due to the 
high gas density and to the reduced size of the system. 

The properties of giant molecular clouds located in the 
galactic centre region can differ significantly from the aver- 
age figures reported above. As an example, we plot in the 
right panel of Fig. [5] the typical time scales for the SgrB2 
cloud. This is a very massive cloud located at 100 pc (pro- 
jected distance) from the galact ic centre. The cloud virial 
mass is M Sg rB2 = 1.9 x 10 6 M Q (|Protheroe et al. 20081 ) and 
a magnetic field at the milliGauss level ha s been measure d 
in the outer envelope of the cloud complex dCrutcher 19991 ). 
The mass distribution can be fit ted with a radial gau ssian 
density profile with a = 2.75 pc (|Protheroe et al. 20081 . and 
references therein) . To compute the curves plotted in Fig. [2] 
(right panel), we assumed a cloud radius of Rs gr B2 = 2<r = 
5.5 pc, which encloses ft; 95% of the total cloud mass. This 
gives an average density of n gas = 1.1 x 10 5 cm~ 3 (roughly 
a factor of 2 below the central density), ft is evident from 
the right panel of Fig. [2] that the SgrB2 cloud is remarkably 
different from a typical giant molecular cloud. In particu- 
lar, the very high values of the magnetic field and of the 
gas density make the energy loss time of CR protons signifi- 
cantly shorter than the propagation time for energies below 
a few hundred GeVs. Moreover, for CR electrons the en- 
ergy loss time is always shorter than the propagation time. 
This means that CR protons with energies up to few hun- 
dred GeVs cannot penetrate the molecular cloud, as they 
do in the cases considered in the left and middle panel of 
Fig. [2] Primary CR electrons cannot penetrate the cloud, 
while secondary CR electrons produced inside the cloud in 
hadronic interactions cannot leave the cloud and radiate all 
their energy close to their production site. These character- 
istics make SgrB2 a very peculiar objects whose modelling 
needs a specific treatment. A detailed study of the CR pene- 
tration in the SgrB2 cloud has been performed by Protheroe 
et al. (2008) , who also computed the synchrotron radio emis- 
sion from secondary electrons. The effects of CRs exclusion 
from giant molecular clouds on their gamma ray emission 
have been discussed in detail by Gabici et al. (2007). In 
the following we will not focus onto any specific object but 
rather investigate the case of the typical molecular clouds as 
the ones considered in the left and middle panels of Fig. [5] 

To conclude the discussion on characteristic time scales, 
it is interesting to note that under certain conditions all the 
energy loss times, for both protons and electrons, scales in 
the same way with gas density, namely as n~^ s . This is ev- 
ident from Eqns. [T5l and [T6l that describe inelastic proton- 
proton scattering and ionization losses for protons respec- 
tively, and from Eqns.[T8]and[l9]that describe ionization and 
Bremsstrahlung losses for electrons respectively. Moreover, 
also for synchrotron losses (Eq. l20p it is possible to derive the 
same scaling with gas density by recalling that observations 
suggest that the magnetic field i n molecular clo uds scales 
as the square root of gas density (|Crutcher 1 999). This has 
the important consequence that, in typical molecular clouds, 
particles with a given energy (protons or electrons) lose their 



energy always through the same channel, independently on 
the gas density. As said above, this conclusion may not hold 
for peculiar obkects such as SgrB2. 



5 NON-THERMAL RADIATION FROM A 
MOLECULAR CLOUD 

Molecula r clouds are now es tabilished gamma ra y 
sources (jAharonian et al. 2006d ; lAharonian et al. 2008al ) . 
Their gamma ray emission is believed to be the result 
of the decay of neutral pions produced during the in- 
elastic collision s of CRs with the dense gas which consti- 
tutes the cloud |lssa and Wolfendale" 19811 ; lAharonian 199ll ; 
iGabici et al. 20071) . During the same interactions, also elec- 
trons and positrons are produced via the decay of charged 
pions. These electrons and positrons produce a broad spec- 
trum of radiation from radio waves to gamma rays due to 
synchrotron emission and non thermal Bremsstrahlung. 

In this section we compute the expected non-thermal 
emission from a molecular cloud located in the proximity of 
a SNR. The emission is the result of CR interactions with 
the dense gas and magnetic field in the cloud and is made 
up of two contributions: a steady state contribution from the 
interactions of background CRs that penetrate the cloud and 
a time dependent contribution from the interactions of CRs 
coming from the nearby SNR. 

We consider a giant molecular cloud of mass M c i = 
1O 5 M0, radius R c i = 20 pc and we assume an uniform den- 
sity of rigas ~ 120 cm -3 . The magnetic field is B c i — 20 /iG. 
The relevant time scales for CR propagation and energy 
losses in such an environment have been plotted in the left 
panel of Fig. [5] In order to show all the different contribu- 
tions to the total non-thermal emission, in Fig. [3] we plot 
the broad band spectrum from the cloud at a time t — 2000 
years after the supernova explosion. The SNR is located 100 
pc away from the molecular cloud and the distance to the 
observer is 1 kpc. The dotted line (curve 3) represents the 
emission from neutral pion decay (from both background 
CRs and CRs from the SNR), the dot-dashed lines repre- 
sent the synchrotron (curve 2) and Bremsstrahlung (curve 
4) emission from background CR electrons that penetrate 
the molecular cloud and the dashed lines represent the syn- 
chrotron (curve 1) and Bremsstrahlung (curve 5) emission 
from secondary electrons produced during inelastic CR in- 
teractions in the dense gas. For the spectrum of background 
CR electrons we use a fit to the measured spectrum at Earth 
(see Kobayashi et al. 2004 for a recent compilation of data). 
Electrons can freely penetrate the cloud except for the high- 
est (> 300 TeV) and lowest (< 100 MeV) part of the energy 
spectrum, where the energy loss time scale is significantly 
shorter than the propagation time (see Fig. [2}. For this en- 
ergies the CR electron flux inside the cloud is suppressed 
and we estimated the suppression by assuming that CR elec- 
trons can penetrate undisturbed the cloud only up to a given 
depth, which is estimated as « CTf ogg and ~ */ Dt^ oss for the 
high and low energy end of the spectrum respectively. This 
approximation is satisfactory for the purposes of this paper, 
given that the suppression becomes crucial only for particles 
which emit negligible non-thermal emission. 

The decay of neutral pions dominates the total emis- 
sion for energies above m 100 MeV. The two peaks in the 
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Figure 3. Broad band spectrum for a molecular cloud of mass 10 5 Mq, radius 20 pc, density ~ 120 cm -3 , magnetic field 20/iG. The 
molecular cloud is at 100 pc from a SNR that exploded 2000 yr ago. The distance of the cloud is 1 kpc. The dotted line shows the 
emission from 7r°— decay (curve 3), the dot-dashed lines represent the synchrotron (curve 2) and Bremsstrahlung (curve 4) emission from 
background CR electrons that penetrate the molecular cloud and the dashed lines the synchrotron (curve 1) and Bremsstrahlung (curve 
5) emission from secondary electrons. 



emission reflects the shape of the underlying CR spectrum, 
which, as illustrated in Fig. [T] is the sum of the steep back- 
ground CR spectrum, which produces the 7r°-bump at a 
photon energy of m,o /2 ~ 70 MeV (in the photon flux F), 
and an hard CR component coming from the SNR that pro- 
duces the bump at higher energies. The flux level at 1 TeV is 
approximatively 5x 10~ 12 erg/cm 2 /s, detectable by currently 
operating Imaging Atmospheric Cherenkov Telescopes, even 
taking into account the quite extended (« 2°) nature of 
the source. It is remarkable that such a cloud would be de- 
tectable even if it were located at the distance of the Galactic 
centre, as can be easily estimated by taking into account that 
the sensitivity of a Cherenkov telescope like H.E.S.S. after 
50 hours of exposure, is w 10" 13 (6» s /0.1°)TeV/cm 2 /s, where 
9 S is the source extension. This means that very massive 
clouds can be used to reveal the presence of enhancements 
of the CR density in different locations throughout the whole 
Galaxy. Similar conclusions can be drawn for the expected 
emission in the GeV range, which is currently probed by the 
AGILE and FERMI satellites. In particular, FERMI, with 
a point source sensitivity of < 10 _9 GeV/cm 2 /s at energies 
above 1 GeV (www-glast.slac.stanford.edu), will be able to 
detect such giant molecular clouds as extended sources if 
they are located within 1-^2 kpc from the Earth, or as point 
sources if they are at larger distances. Such a use of molec- 
ular cloud as CR baro meters has been discus sed in several 
papers for b oth GeV illssa and Wolfendale 198lh and TeV 
gamma rays (| Aharonian 1991 ). Here we demonstrated that 
SNRs can provide enhancements in the CR density that can 
generate gamma ray fluxes well within the capabilities of 
currently operating instruments. 

The spectral shape in the gamma ray energy range de- 
serves further discussion. For the situation considered in 
Fig. [3] the GeV gamma rays are the result of the decay of 
neutral pions produced by background CRs that penetrate 
the cloud. Thus, the gamma ray spectrum above GeV ener- 



gies simply mimics the underlying CR spectrum, which is a 
steep power law of the form ~ E~ 2,75 . On the other hand, 
the neutral pion decay spectrum at TeV energies is, in this 
case, dominated by the contribution from CRs coming from 
the nearby SNR. After 2000 years from the supernova explo- 
sion, only CRs with energies above several tens of TeV had 
enough time to leave the SNR and reach the cloud and thus 
the CR spectrum at the cloud exhibits an abrupt low energy 
cutoff at that energy, that we call here (see Fig. [TJ . As 
a consequence, the gamma ray spectrum is expected to be 
peaked at an energy ~ E^^/IO of several TeV. The slope 
of the gamma ray spectrum below the peak is determined 
by the physics of the interaction only, and not by the shape 
of the underlying CR spectrum, and is r oughly of the form 
E 2 dN/dE oc E 2 (e.g. lKelner et al. 20061 ). Thus, a loose as- 
sociation between a SNR and a massive molecular clouds 
as the one studied here, is expected to be characterized, at 
least at some stage of the SNR evolution, by a very peculiar 
spectrum which is steep at low (GeV) energies and hard at 
high (TeV) energies. 

The possibility of detecting sources with such a distinct 
spectrum is also relevant for the issue of identifying GeV and 
TeV unidentified sources. One of the criteria suggested to 
support an association between a GeV and TeV source is, be- 
side the positional coincidence, the spectral compatibility. In 
a recent paper, Funk et al. (2008), investigated the spectral 
compatibility of EGRET and H.E.S.S. unidentified sources 
located in the inner Galactic region. For sources showing 
positional coincidence, they found generally a good spectral 
compatibility, but due to the small number of sources, it 
was not possible to claim an association between the two 
populations of sources at a statistically significant level. For 
sources detected by only one of the two instruments, Funk 
et al. (2008) investigated the consequences of extrapolat- 
ing the measured spectrumat higher or lower energies. They 
considered only single power laws extending over the whole 
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Figure 4. Time dependence of the broad band spectrum for the same giant molecular cloud considered in Fig. [3] The solid, dotted, 
short-dashed and long-dashed lines refers to the emission at a time 500, 2000, 8000 and 32000 years after the explosion. The dot-dashed 
line represent the spectrum of a molecular cloud with no SNR in its proximity. 



(GeV-TeV) energy range and discussed the implications of 
adding a high or low energy cutoff. The scenario presented 
here adds the new interesting possibility of a single source 
showing a dramatically different behavior at GeV and TeV 
energies, namely a spectrum which shows a significant hard- 
ening at higher energies. 

The evolution with time of the emission from the cloud 
is shown in Fig. [4] where the solid, dotted, short-dashed and 
long-dashed lines show the spectrum at a time equal to 500, 
2000, 8000 and 32000 years after the supernova explosion 
respectively. For comparison, the emission from a molecular 
cloud, with no SNR located in its proximity, is plotted as 
a dot-dashed line. In this case, only background CRs that 
penetrate the molecular cloud contribute to the emission. 

It is evident from Fig. [4] that the radio (A > 0.1mm) 
and the soft gamma ray (~ f MeV -rl GeV) emission from 
the cloud is constant in time. This reflects the fact that 
such emission is produced by background CRs that enter the 
cloud. The emission in the other energy bands is variable in 
time, being produced by the CRs coming from the SNR. The 
flux of this latter component changes with time as indicated 
in Fig. [T] The two most prominent features in the variable 
emission are two peaks, in X- and gamma-rays respectively. 

The peak observable at gamma ray energies is the re- 
sult of the decay of neutral pions produced when CRs of 
different energies coming from the SNR reach the cloud at 
different times. The peak moves at lower and lower energies 
with time, reflecting the fact that CRs with lower and lower 
energies progressively reach the cloud as time flows. At early 
times, the emission can extend up to ~ 100 TeV (solid line), 
revealing the presence of PeV CRs and thus indirectly the 
fact that the nearby SNR i s acting as a CR pevatron (see 
iGabici and Aharonian 20071 . for a discussion of this issue). 
Moreover, the gamma ray emission in the TeV range is en- 
hanced with respect to the one expected from an isolated 
molecular cloud (dashed-dotted line in Fig. [4} for a period 
of several 1 4 yr. This is much longer than the period dur- 
ing which SNRs are effectively accelerating the multi-TeV 
CRs responsible for the TeV emission, which lasts few thou- 
sands years. This is because the duration of the gamma ray 



emission from the cloud is determined by the time of prop- 
agation of CRs from the SNR to the cloud and not by the 
much shorter CR confinement time in the SNR. Therefore, 
the gamma ray emission from the cloud lasts much longer 
than the emission fro m the SNR, making the dete ction of 
clouds more probable (jGabici and Aharonian 2007h . 

The peak in the X-ray spectrum is due to synchrotron 
emission from secondary electrons produced in CR interac- 
tions in the cloud. The peak is moving to lower energies 
with time but, unlike the gamma ray peak, it is also becom- 
ing less and less pronounced. This fact can be understood 
with the help of Fig. (left panel) . After 500 yr from the 
supernova explosion (solid line in Fig. [3]), PeV CRs from 
the SNR reach the cloud and produce there secondary elec- 
trons with energy in the « 100 TeV range. For these elec- 
trons, the synchrotron cooling time is comparable with the 
escaping time from the cloud. Thus, they release a consid- 
erable fraction of their energy in form of X-ray synchrotron 
photons before leaving the cloud. As time passes, lower en- 
ergies CRs reach the cloud and secondary electrons with 
lower energies are produced. For these electrons the cool- 
ing time becomes progressively longer than the escape time 
and this explain the suppression of the synchrotron emis- 
sion. The X-ray synchrotron emission is weaker than the 
TeV emission for any time and for times > 2000 yr the 
ratio between TeV and keV emission can reach extreme val- 
ues of a few tens or more. These values are observed from 
some of the unidentified TeV sources such as H ESS J1616- 
508 (|Matsumoto et al. 20071 : iBamba et al. 20071 ) and more 
in general unidentified TeV sources are characterized by the 
absence of any clear counterpart at other wavelengths. Due 
to this peculiar spectral properties, such sources have been 
labeled as dark, since they seemed to emit gamma rays only. 
However, in the scenario presented in this paper, spectra 
showing a high TeV/kev flux ratio can be produced very 
naturally if a cloud is illuminated by CRs coming from a 
nearby SNR. This suggestion is also supported by the fact 
that most of the unidentified TeV sources are spatially ex- 
tended, as molecular clouds are expected to be. 

In the hard X-ray /soft gamma ray region of the spec- 
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trum (from tens of keVs to hundreds of MeVs), partially 
covered by the INTEGRAL satellite, extremely hard spec- 
tra (dN/dE oc E -1 ) may result due to the Bremsstrahlung 
emission from primary and secondary CR electrons. 

Finally, the radio emission from the cloud is the re- 
sult of the synchrotron emission of background CR elec- 
trons that penetrate the cloud. The contribution from sec- 
ondary electrons is subdominant, due to the fact that the 
~ GeV secondary electrons that might emit synchrotron 
radio waves escape from the cloud before losing energy 
(see Fig. [2}. Moreover, in the GeV energy range, ioniza- 
tion and Bremsstrahlung losses dominate, and this would 
further reduce the expected synchrotron emission. Recently, 
Protheroe et al. (2008) estimated the expected synchrotron 
radio emission from secondary electrons produced by CR 
interactions in the Sgr B2 giant molecular cloud. They did 
not consider the contribution from background CR elec- 
trons and neglected the effects of diffusive transport of sec- 
ondary electrons (namely, they assumed no penetration of 
electrons from outside the cloud and instantaneous cooling 
of secondary electrons produced inside the cloud) . These as- 
sumptions are valid for the very strong magnetic field of 
the order of one milliGauss measured for SgrB2. Such field 
is much stronger than the one assumed here. Thus, due to 
these somewhat extreme assumptions, possibly justified in 
the particular case of Sgr B2 cloud, a direct comparison be- 
tween their findings and the results presented in this paper 
is not straightforward. Protheroe et al. (2008) also noticed 
that the radio emission from a cloud at frequencies above 
~ 1 GHz can be dominated by free-free thermal emission. 
Being focused on the non-thermal emission from clouds, we 
did not attempt here to model their thermal emision. How- 
ever, our findings can still be tested and constrained by ob- 
servations in the GHz range. This can be done by requiring 
the predicted synchrotron emission not to overproduce the 
observed free-free thermal emission. A similar approach has 
been adopted by Jones et al. (2008). 

5.1 The spectral shape at GeV-TeV energies 

As noticed in the previous section, concave gamma ray spec- 
tra may be produced in a molecular cloud located in proxim- 
ity of a SNR, as the result of the decay of neutral pions pro- 
duced in CR interactions. Such concavity reflects the shape 
of the underlying CR spectrum, which consist of the super- 
position of two components: the galactic CR background, 
characterized by a steep spectrum, and the CRs coming from 
the nearby SNR, which exhibit a hard spectrum. With this 
respect, the distance between the SNR and the molecular 
cloud d c i plays a crucial role. This is because, the larger the 
distance between the SNR and the cloud, the lower the level 
of the CR flux coming from the SNR. Moreover, also the 
time evolution of the emission from a cloud changes with 
d c i since the time it takes a particle with given energy to 
cover such a distance scales as t ~ d^i/D, where D is the 
diffusion coefficient. 

In Fig.[S]the total gamma ray spectrum from a molecu- 
lar cloud is shown as a function of the distance between the 
SNR and the cloud. The cloud mass is 10 5 Mq and the dis- 
tance from the SNR is 50, 100 and 200 pc for the left, central 
and right panel, respectively. Similarly to Fig. [4] the solid, 
dotted, short-dashed and long-dashed lines refer to the emis- 



sion for 500, 2000, 8000 and 32000 years after the supernova 
explosion. It is evident from Fig. [5] that a great variety of 
gamma ray spectra can be produced. In almost the entirety 
of the cases considered, the gamma ray emission is charac- 
terized by the presence of two pronounced peaks. The low 
energy peak, located in the GeV domain is steady in time 
and it is the result of the decay of neutral pions produced 
in hadronic interactions of background CRs in the dense in- 
tercloud gas. The high energy peak is the result of hadronic 
interactions of CRs coming from the nearby SNR, and thus 
it is moving in time to lower and lower energies (see previous 
section for a discussion) . Both the relative intensity and po- 
sition of the two peaks depend on the distance between the 
SNR and the cloud. Interestingly, the GeV emission from the 
cloud is affected by the presence of the nearby SNR only at 
late times after the explosion and only if the distance from 
the SNR is comparable or smaller than ps 50 pc (see Fig. [5] 
left panel) . In all the other cases the GeV emission is always 
the result of the interactions of background CRs and thus, 
at least in this case, observations of molecular clouds in the 
GeV gamma ray domain cannot be used to infer the pres- 
ence of a CR accelerator located at a distance greater that 
~ 50 pc from the cloud. 

Similar concave or "V-shaped" spectra have been re- 
cently obtained by Rodriguez Marrero et al. (2008) in a dif- 
ferent context in which two molecular clouds are assumed 
to be located in the proximity of a CR accelerator. If the 
two clouds happen to be located at different distances from 
the CR accelerator but within an angular separation smaller 
than the FERMI angular resolution, then they would appear 
as a single GeV source, and the superposition of their emis- 
sion might result in concave spectra. However, Rodriguez 
Marrero et al. (2008) did not include in their calculations 
the contribution to the emission coming from the ubiquitous 
galactic CR background, which in most cases dominates over 
the contribution of CRs from the nearby SNR, at least for 
what concerns the GeV emission from the cloud (see e.g. 
Fig. O . Moreover, they also considered very short distances 
between the cloud and the CR accelerator (down to 5 pc), 
but in this case an accurate modeling of the CR accelerator 
itself is needed, especially for what concerns its own gamma 
ray emission that might add up to the total emission. In ad- 
dition to that, distances as short as ~ 5 -f- 10 pc are, under 
many circumstances, smaller than both the source size and 
the radius of the molecular cloud itself, and this changes 
significantly the problem since the interactions between the 
accelerator and the MC are likely to play an important role. 

The prediction of V-shaped spectra that we make in this 
paper is more general than the one by Rodriguez Marrero 
et al. (2008), since it represents an intrinsic feature of a 
single molecular cloud which is located close to a CR source. 
The V-shaped gamma ray spectrum reflects the shape of 
the underlying CR spectrum which is the superposition of 
the steep spectrum of the background CRs and the hard 
spectrum of CRs coming from the nearby SNR. 

5.2 The role of the magnetic field 

The value of the magnetic field in the cloud is an important 
parameter since it regulates the synchrotron energy losses of 
high energy electrons, and the diffusive escape time of rela- 
tivistic particles from the cloud. Observations suggest that 
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Figure 5. Total gamma ray emission from a molecular cloud of mass 10 Mq located at a distance of 1 kpc. The distance between 
the molecular cloud and the SNR is 50, 100 and 200 pc for left, centre and right panel, respectively. As in Fig. [4] the solid, dotted, 
short-dashed and long— dashed lines refers to the emission at a time 500, 2000, 8000 and 32000 years after the explosion. 



the value of the magnetic field in a molecular cloud scales 
with the square root of the gas density, reaching very high 
values (1 m G or more) in t he very dense sub-parsec scale 
cloud cores |Crutcher 19991 ). However, dense cores consti- 
tute a very small fraction of the total volume of molecular 
clouds. Since here we are interested in calculating the emis- 
sion from the whole cloud, the relevant parameter is the 
volume averaged value of the cloud magnetic field. It seems 
reasonable to assume that the total magnetic energy in the 
cloud Wb = (B 2 /8n) x V, where V is the cloud volume, does 
not exceed the total gravitational energy of the cloud Wg = 
SGM 2 /5R ct . This leads to a maximum value of the aver- 
age magnetic field of: B sC 30 (M/10 5 M Q )(i? ci /20pc)- 2 /iG- 
This is in general agreement with observations, from which 
a value of « 10 /j,G can be inferred for typical cloud den- 
sities of « 100 cm -3 . However, since the dispersion aro und 
this mean value is considerable (see e.g. ICrutcher 1999T ). it 
is worth to investigate how the non thermal emission from a 
molecular cloud depends on the actual value of the magnetic 
field. 

Fig. |5]shows the broad band spectrum from a molecular 
cloud of mass 1O 5 M0, radius 20 pc, density ~ 120 cm -3 . The 
distance between the cloud and the SNR is 50 pc. In the top 
and bottom panels is plotted the cloud emission after 500 
and 2000 years from the supernova explosion, respectively. 
The solid line refers to a value of the magnetic field of 30 /iG, 
while the dashed line represents the emission for a smaller 
value of the field equal to 10 /iG. The cloud emission from 
radio frequencies up to the hard X-ray band strongly de- 
pends on the value of the magnetic field, while the gamma 
ray emission is unaffected, being the result of hadronic inter- 
actions of CR protons. The strong dependence of the radio 
and X-ray emission on the magnetic fields is evident, and 
this demonstrate that this effect has to be taken into ac- 
count in multi wavelength studies of molecular clouds. 



6 CONCLUSIONS 

The non thermal emission from a molecular cloud located 
in proximity of a SNR is the result of the interactions of 
CRs that penetrate the cloud. Both CRs in the galactic 



-10 




Log(E) [keV] 

Figure 6. Broad band emission from a molecular cloud of mass 
10 5 Mq, radius 20 pc, density ~ 120 cm -3 . The distance between 
the cloud and the SNR is 50 pc. Top and bottom panels refer 
to 500 and 2000 years after the supernova explosion, respectively. 
The magnetic field is 30 fiG (solid lines) and 10 fiG (dashed lines). 

background and runaway CRs from the SNR contribute to 
the emission. In this paper, we calculated the expected non 
thermal emission from the molecular cloud as a function of 
several parameters, such as the time after the supernova ex- 
plosion, the distance between SNR and cloud, and the cloud 
magnetic field. In all the cases, the gamma ray emission from 
the cloud is by far exceeding the energy output in other en- 
ergy bands. 

The gamma ray emission from the cloud consists of two 
distinct components. The first one is due to the interactions 
of CRs from the galactic background, which are character- 
ized by a steep spectrum. This component is steady in time 
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and peaks in the GeV energy region. The second component 
is the result of the interactions of runaway CRs that escape 
the SNR and diffusively reach the cloud. The spectrum of 
these runaway CRs is hard and variable in time and both 
this characteristics are reflected in the related gamma ray 
emission. In particular, this component is producing a sec- 
ond peak in the gamma ray spectrum which moves to lower 
and lower energies, reflecting the fact that high energy CRs 
are released from the SNR earlier and diffuse faster than 
low energy ones, and as a consequence they reach the cloud 
earlier. The superposition of these two components of the 
gamma ray emission may produce concave spectra, which 
are steep at low (~ GeV) energies and hard at high (~ 
TeV) energies. This peculiar spectra might be revealed by 
means of joint observations of FERMI and ground based 
Cherenkov telescopes. The detection of such emission would 
provide a strong hint for the presence of a CR accelerator 
in the neighbourhood of the molecular cloud. 

It has been suggested that some of the dark TeV 
unidentified sources, with no obvious counterparts at 
any other wavelength, might be indeed associated with 
molecular clouds illuminate d by CRs from nearby ac- 
celerators such as SNRs l|Aharonian and Atovan 19961 ; 
iGabici and Aharonian 20071 ). We showed here that, in this 
scenario, the "darkness" of a source, often denned as the 
radio Rrev/x between the TeV and X-ray observed flux, 
depends on several parameters such as the distance between 
the SNR and the cloud, the time after the supernova ex- 
plosion and the value of the magnetic field in the molecular 
cloud. However, in most cases it is very natural to obtain 
very high values for RreV/Xi compatible with the estimates 
on RteV/ x claimed for some of the unidentified T eV sources 
(see e.g. iBamba et al. 20071 : iBamba et al. 20081 ). This fur- 
ther supports an association between unidentified TeV 
sources and molecular clouds and suggests that the unknown 
nature of these objects might be soon revealed. 
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